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ABSTRACT 


The  primary  goals  of  this  thesis  included  developing  a  computer-integrated  LDV 
system,  utilizing  the  LDV  system  to  measure  the  mean  and  fluctuating  quantities  of  the 
flow  field  of  an  airfoil-spoiler  system,  and  comparing  the  results  of  the  experiment  with 
previous  hot-wire  anemomctry  results. 

The  experiment  was  conducted  in  the  NASA  Ames  1 1  cm  x  25  cm  indraft 
tunnel.  A  two-component  LDV  system  was  developed  and  set  up,  and  measurements 
were  obtained  along  the  upper  surface  of  the  airfoil  (ahead  of  and  behind  the  spoiler) 
as  well  as  in  its  wake,  close  to  the  trailing  edge.  Mean  velocity  profiles  and  turbulent 
stresses  were  calculated  and  comparisons  were  made  with  previous  hot-wire 
experiments. 

It  was  shown  that  the  hot-wire  technique  is  not  as  accurate  as  the  LDV  in  highly 
turbulent  regions  of  the  flow  field.  The  hot-wire  results  are  also  not  reliable  in  the  near 
wake  region  due  to  the  inherent  reverse  flow.  In  other,  less  turbulent,  regions  of  the 
(low  field  where  there  is  no  reverse  flow,  good  agreement  between  previous  hot-wire 
data  and  present  experimental  results  was  found. 
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I.  INTRODUCTION 


A.  SPOILERS 

A  spoiler  is  a  plate  on  the  upper  surface  of  a  wing  which  can  be  deflected  to  spoil 
the  flow  and  to  cause  early  separation.  Consequently,  lift  decreases  and  drag  increases. 
Spoilers  are  used  as  effective  aerodynamic  control  surfaces  in  transport  aircrafts.  They 
can  be  symmetrically  deflected  as  speed  brakes  or  as  lift  dumpers  at  touch  down. 
When  spoilers  are  deflected  asymmetrically,  they  can  provide  substantial  amounts  of 
roll  control  in  certain  conditions  of  flight. 

The  flow  field  of  an  airfoil  with  a  deflected  spoiler  is  complex:  including  flow 
separation,  reattachmcnt  and  vortex  shedding  [Ref.  l:p.  1], 

The  basic  steady  state  structure  of  the  flow  field  of  an  airfoil  with  a  deflected 
spoiler  has  been  measured  and  hypothesized  by  the  Boeing  Research  Group  [Ref.  2). 
Further  work  into  how  tne  unsteady  wake  relates  to  the  airfoil  geometry  and  to  other 
parameters  was  undertaken  by  Lee.  ct  al.  In  their  experiments,  Lee,  et  al.,  used 
hot-wire  anemometers  and  pitot  tubes  to  systematically  measure  the  mean  and 
fluctuating  quantities  on  the  airfoil  surface,  as  well  as  in  the  wake,  in  order  to  better 
understand  the  structure  of  this  How  field.  [Ref.  1| 

B.  LASER  DOPPLER  VELOCIMETER  SYSTEMS 

Advanced  research  in  experimental  fluid  dynamics  requires  a  familiarity  with 
sophisticated  measurement  techniques.  Optical  methods  and  in  particular,  the  Laser 
Doppler  Vcloeimeter  (LDV)  are  now  recognized  as  the  most  reliable  means  for 
performing  measurements  in  complex  turbulent  flows  [Ref.  3]. 

Since  its  inception  in  1964,  the  LDV  technique  has  been  revised,  developed,  and 
tested  by  many  researchers  and  its  capabilities  arc  reasonably  well-understood,  flic 
method  has  the  advantage  of  performing  instantaneous,  non-intrusive  velocity 
measurements  of  small  particles  suspended  in  the  flow.  These  measurements  arc 
relatively  independent  of  fluid  properties  and  the  method  docs  not  require  frequent 
calibration.  Furthermore,  its  ability  to  measure  the  flow  without  disturbing  it  is 
especially  attractive.  However,  because  of  the  range  of  disciplines  involved  including 
optics,  electronics,  light  scattering,  and  signal  processing,  the  implementation  of  this 
method  can  be  quite  onerous.  Serious  measurement  errors  can  occur  if  the 
measurement  techniques  are  not  applied  correctly. 


C.  THESIS  GOALS 

Considering  this  previous  work,  the  goals  of  this  thesis  experiment  were  threefold. 
First,  to  develop  and  set  up  a  two-component  LDV  system  with  complete 
computer-experiment  integration  employed.  Second,  to  utilize  this  LDV  system  to 
measure  the  mean  and  fluctuating  velocities  of  an  airfoil  with  a  deflected  spoiler  flap, 
especially  in  the  reverse  flow  region  and  the  near  wake  where  the  fluctuations  are  large. 
Third,  to  compare  the  LDV  data  obtained  with  the  previous  hot-wire  anemometry 


The  experiments  were  conducted  in  the  NASA  Ames  1 1  cm  x  25  cm  indraft 
wind  tunnel.  Figure  2.1  shows  the  test  section  of  this  tunnel.  The  test  section  is 
equipped  with  slotted  walls  and  a  number  of  plenum  chambers  above  and  below  the 
centerline.  These  features  were  used  in  previous  experiments  to  give  the  tunnel 
adaptive  walls.  The  adaptive  wall  concept  was  not  employed  in  this  experiment, 
however  the  slotted  walls  did  provide  a  means  of  reducing  wall  interference.  The  wind 
tunnel  also  is  equipped  with  interchangeable  glass  and  plcxiglas  test  section  windows. 
For  these  experiments,  the  glass  windows  were  used  due  to  their  higher  refractive 
qualities  which  enhance  the  LDV  signal*. 

The  basic  model  used  was  a  Boeing  advanced  transport  research  airfoil  section 
with  a  6-inch  chord  and  maximum  thickness  ratio  of  11.3%  as  shown  in  Figure  2.2. 
This  airfoil  has  a  spoiler  of  15.5%  chord  length  hinged  at  73%  chord.  The  coordinates 
of  the  airfoil  are  given  in  Table  I.  This  airfoil  was  selected  because  of  the  availability 
of  hot-wire  anemometer  and  pitot  tube  experimental  data  for  comparison  with  the 
two-component  LDV  results. 

B.  INSTRUMENTATION 

1.  Pressure  Instrumentation 

For  this  experiment,  pressure  information  was  obtained  by  taking  advantage 
of  a  bank  of  six,  24-port  Scanivalves  tapped  into  the  tunnel  during  a  previous 
experiment.  The  six  pressure  transducers  were  of  the  differential  type  with  a  range  of 
±5  psig.  One  side  of  each  transducer  was  connected  with  the  measured  pressure 
through  plastic  tubing  while  the  other  side  of  each  transducer  was  vented  to 
atmosphere. 

In  order  to  calibrate  the  transducers,  one  port  of  each  Scanivalve  was  tapped 
to  a  calibration  pressure  generated  by  a  small  electric  vacuum  pump.  Another  port  of 
each  Scanivalve  was  tapped  to  the  atmosphere.  The  calibration  pressure  and 
atmospheric  pressure  were  each  measured  by  alternate  methods.  These  two  known 
pressures  and  the  values  given  for  them  by  each  transducer  allowed  the  Scanivalves  to 
be  calibrated. 


vvs  s  s  VV 


Figure  2.2  Boeing  advanced  transport  research  airfoil  section. 

2.  Laser  Doppler  Velocimeter  Instrumentation 

Figure  2.3  shows  the  two-component,  two  color  LDV  arrangement  used  at 
NASA  Ames.  The  laser  is  a  four  Watt  Argon  ion  unit  (Spectra  Physics  Model  64). 
This  laser  is  configured  with  an  end  cavity  multiline  reflector  which  allows  it  to  operate 
in  the  multiline  mode.  Translation  stages  in  three  axes  permit  positioning  of  the  probe 
volume  anywhere  in  the  test  section. 

In  the  multiline  mode,  the  resultant  beam  is  comprised  of  all  wavelengths 
inherent  to  the  Argon  ion  unit.  This  beam  exits  the  laser  in  a  fairly  collimated  state, 
hence  an  external  collimator  is  not  required.  The  beam  first  passes  through  a 
half-wave  plate  (polarization  rotator)  where  it  becomes  horizontally  polarized.  A  set  of 
Brewster  angle  prisms  is  then  used  to  separate  the  primary  beam  into  its  various 
wavelengths.  To  obtain  maximum  light  transmission  through  the  prisms,  the  beam 
polarization  must  be  parallel  to  the  plane  defined  by  the  incident  and  refracted  beams 
[Ref.  4],  This  was  the  reason  that  the  initially  vertically-polarized  laser  beam  was 
passed  through  a  polarization  rotator  prior  to  entering  the  dispersion  prism.  When  the 
two  prism  method  is  used,  the  prisms  also  turn  the  beams  which,  with  a  little  planning, 
may  be  used  to  advantage.  At  some  distance  from  the  prisms,  the  beams  have 


TABLE  I 

AIRFOIL  COORDINATES  (INCHES) 

X 

Y 

upper 

Y 

1  lower 

0.0 

0.0 

0.0 

.0030 

.0222 

-.0108 

.0060 

.0300 

-.0162 

.0150 

.0468 

-.0258 

.0300 

.0672 

-.0348 

.0900 

.1242 

-.0588 

.1200 

.1452 

-.0672 

.  1 500 

.1626 

-.0750 

.2250 

.2004 

-.0912 

.3000 

.2304 

-.1050 

.4500 

.2790 

-.1296 

.6000 

.3126 

-.1524 

.7500 

.3396 

-.1728 

.9000 

.3594 

-.1920 

1.2000 

.3858 

-.2250 

1.5000 

.3990 

-.2502 

l.sooo 

.4044 

-.2670 

2.1000 

.4038 

-.274S 

2.4000 

.4002 

-.2742 

2.7000 

.3924 

-.2658 

3.0000 

.3816 

-.2502 

3.3000 

.367S 

-.2298 

3.6000 

.3510 

-.2064 

3.9000 

.3288 

-.1818 

4.2000 

.3018 

-.1560 

4.5000 

.2700 

-.1308 

4.8000 

.2304 

-.1044 

5.1000 

.1842 

-.0792 

5.4000 

.1308 

-.0540 

5.7000 

.0696 

-.0282 

6.0000 

.0024 

-.0024 

sufficient  separation  so  mirrors  can  be  used  to  direct  the  individual  beams  to  the  next 
optical  components.  The  two  stronger  lines,  green  (514.5  nm)  and  blue  (488  nm),  are 
guided  by  means  of  turning  mirrors  to  Bragg  cells.  The  weaker,  unused  lines  arc 
"trapped"  (blocked  by  nonreflective  plates)  to  keep  them  from  interfering  with  the 
system.  Prior  to  reaching  its  Bragg  cell,  the  green  beam  is  passed  through  another 
half-wave  plate  in  order  to  vertically  polarize  it. 

A  Bragg  cell  is  an  acousto-optic  device.  Both  the  blue  and  green  beams  are 
passed  through  Bragg  cells  which  split  each  into  two  equal-intensity,  slightly  diverging 
beams  while  inducing  a  40  MIlz  modulation  in  one  of  the  beams  in  each  pair.  The 
dilference  is  that  the  green  beam  is  split  in  the  vertical  plane  while  the  blue  beam  is 
split  in  the  horizontal  plane  (Figure  2.3). 


In  each  case  (blue  and  green),  the  two  diverging  beams  are  then  made  parallel 
by  passing  each  pair  through  its  own  alignment  prism.  From  there  on,  turning  mirrors 
guide  the  beams  through  the  translation  stages.  They  are  finally  focused  by  means  of  a 
condensing  lens  at  the  desired  point  in  the  test  section  of  the  wind  tunnel. 

The  four  beams  must  all  focus  to  a  common  point  in  the  flow  and  be 
approximately  coaxial.  Using  the  steering  mirrors  for  the  appropriate  beam  pair,  the 
focused  beams  are  adjusted  to  visually  overlap.  Following  this,  various  techniques 
exist  for  making  fine  adjustments  to  insure  the  beams  all  cross  at  a  common  point. 

An  off-axis,  forward  scatter  configuration  has  been  selected  for  the  receiving 
optics  used  at  NASA  Ames.  This  gives  flexibility  to  the  system  so  it  can  be  adjusted  to 
different  requirements;  if  maximum  signal  strength  is  required,  angular  deviation  from 
the  beam  axis  of  symmetry  is  minimized.  If  focusing  on  the  smallest  part  of  the  probe 
volume  is  of  interest,  the  angular  deviation  is  made  as  large  as  the  physical  constraints 
will  allow. 

The  first  lens  in  the  receiving  optics  is  the  aperture  through  which  the 
scattered  light  is  collected  (Figure  2.4).  The  second  lens  collimates  the  beam.  A 
collimated  beam  is  required  by  the  polarization  beamsplitter  for  maximum  efficiency. 
The  polarization  beamsplitter  allows  the  vertically  polarized  green  beam  to  pass 
straight  through  while  the  horizontally  polarized  blue  beam  is  reflected  upward.  The 
use  of  polarization  rather  than  frequency  to  separate  the  two  wavelengths  substantially 
increases  the  efficiency  of  the  system.  Each  beam  then  passes  through  a  line  filter 
centered  on  the  respective  wavelength  (blue  or  green)  in  order  to  reduce  both  the 
ambient  light  and  the  cross-talk  due  to  impure  polarization.  Each  line  filter  is  followed 
by  a  lens  which  serves  to  focus  the  light  onto  the  pinhole  aperture  of  the  respective 
photomultiplier.  The  photomultipliers  convert  this  light  energy  into  electronic  signals 
which  are  then  passed  on  to  the  signal  processors. 

Preamplifiers  are  located  in  the  photodetector  housings  to  further  increase  the 
signal  amplitude  and  drive  the  signal  over  relatively  long  cables.  The  preamplifier  is 
kept  close  to  the  detectors  in  order  to  increase  the  signal  level  before  additional  noise 
enters  the  signal  cables.  It  is  always  good  practice  to  amplify  the  signal  as  much  as 
possible  with  the  photomultiplier  tube  before  going  to  the  preamplifier  which  inevitably 
produces  electronic  noise. 

The  signal  processors  used  in  this  experiment  were  of  the  counter  type  and 
produced  by  Macrodync.  These  processors  have  good  resolution  and  superior  noise 
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Figure  2.4  Computer-experiment  integration. 
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rejection  circuitry.  The  Macrodyne  unit  is  a  high-speed  counter  with  a  1  GHz  clock 
designed  to  operate  in  a  variety  of  experimental  situations.  It  has  been  applied  under 
single  particle,  multi-particle  and  multi-velocity  conditions.  In  addition  to  the  basic  5/8 
(or  10/16)  signal  validation  criteria  and  oversize  particle  rejection,  a  sophisticated 
multi-level  sequence  validation  has  been  incorporated.  This  allows  the  experimenter  to 
select  the  criteria  which  will  best  yield  data  for  a  given  signal-to-noise  ratio,  signal 
level,  seeding  level,  Doppler  frequency,  and  fluctuation,  indicative  of  the  test  situation. 

When  the  particle  seeding  rate  is  low,  as  was  the  case  in  this  experiment,  the 
LDV  signals  become  intermittent  bursts.  Under  these  conditions,  the  counter 
processor  is  the  best  choice  as  a  means  for  obtaining  the  Doppler  difference  frequency. 
If  the  seeding  is  high,  the  signals  will  be  almost  continuous  with  random  amplitude  and 
phase  fluctuations.  The  counter  will  still  function  but  cannot  be  expected  to  produce 
accurate  results  to  better  than  ±  1/2  cycle. 

The  counter  processor  determines  the  Doppler  difference  frequency  in  the  time 
domain.  Timing  over  a  selected  number  of  cycles  in  the  burst  signal  is  accomplished 
with  a  crystal  controlled  "clock"  oscillator.  The  frequency  of  the  oscillator  must  be 
much  higher  than  the  Doppler  difference  frequency  to  achieve  good  resolution.  The 
clock  count  is  initiated  when  the  high-pass  filtered  signal  exceeds  a  set  threshold  level. 
A  preselected  number  of  cycles  (8  or  16)  are  processed  and  used  to  stop  the  counter. 
Greater  accuracy  may  be  obtained  by  processing  more  cycles.  Some  processors,  (c.g., 
Acrometrics)  count  all  of  the  cycles  in  the  burst  and  average  the  result  over  the  number 
of  cycles  counted.  This  capability  is  especially  important  when  the  seed  particles  arc 
polydisperse  and  produce  a  range  of  signal  amplitudes.  In  such  cases,  there  is  the 
possibility  of  completing  the  preset  number  of  cycles  on  the  early  part  of  the  signal 
which  has  a  lower  signal-to-noise  ratio  and  slew  rate.  There  is  also  the  possibility  of 
reading  signals  from  larger  particles  more  than  once. 

Although  the  importance  of  attaining  good  quality  signals  cannot  be 
overemphasized,  in  many  practical  cases  the  signals  will  be  noisy.  Therefore,  each 
burst  must  pass  certain  criteria  before  it  is  accepted.  The  so-called  periodicity  check  is 
used  to  reject:  1)  signals  with  insufficient  cycles,  2)  signals  produced  by  two  or  more 
particles  in  the  measuring  volume  at  a  time,  and  3)  noise.  Two  counters  are  used  to 
make  period  measurements  over  a  different  number  of  cycles,  usually  5  and  8.  The 
ratio  of  the  two  measurements  is  then  compared  to  five-eighths  and  if  they  agree  to  a 
specified  percentage,  the  signal  is  validated.  The  problem  with  this  technique  is  that 
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the  entire  signal  may  be  rejected  because  of  a  noisy  cycle  early  in  the  burst.  However, 
when  the  method  is  coupled  with  the  three-level  validation  circuitry,  it  can  be  effective 
in  rejecting  spurious  noise. 

The  three-level  validation  requires  that  a  signal,  Figure  2.5,  exceed  a  positive 
threshold  level,  pass  through  zero,  and  exceed  a  negative  threshold  level  before  the  next 
zero  crossing  is  counted  as  a  cycle.  Such  logic  has  proven  to  be  very'  effective  in 
rejecting  noise.  If  the  conditions  are  not  met  on  any  cycle,  the  counter  resets  and 
starts  again  almost  immediately.  Thus,  noisy  cycles  on  the  start  of  a  burst  need  not 
cause  the  loss  of  an  otherwise  good  signal. 

These  amplified  electronic  signals  are  sent  to  the  signal  processor  and 
ultimately  to  the  data  processor  for  conversion  into  velocities. 

Ideally,  an  LDV  needs  no  calibration  since  it  provides  a  direct  measure  of  the 
particle  velocity.  The  fringe  spacing,  6,  can  be  determined  using  the  following  formula 
where  X  is  the  laser  wavelength  and  0  is  the  angle  of  the  intersecting  beams: 

6  =  X/(2  x  sin(0/2)) 

Knowing  the  fringe  spacing  and  the  time  required  for  a  given  number  of  crossings, 
particle  velocity  can  be  determined.  In  practice  this  method  can  be  inaccurate  due  to 
misalignment  of  the  optics.  The  beam  intersection  angle,  0,  is  also  very  difficult  to 
compute  accurately.  Therefore,  it  is  wise  to  test  the  entire  system  by  measuring  the 
light  scattered  from  an  object  moving  at  a  velocity  that  is  known  a  priori.  A  rotating 
disk  with  a  mark  at  a  known  radius  is  usually  the  most  convenient  means  of  doing  this. 
However,  another  method  (the  method  used  in  this  experiment)  is  to  measure  the  free 
stream  velocity  using  the  LDV  and  compare  this  value  with  that  computed  from  the 
Scanivalve  data. 

3.  Flow  Seeding  Requirements 

An  inconvenient  characteristic  of  the  laser  Doppler  velocimctcr  technique  is 
the  need  for  seed  particles.  The  LDV  actually  measures  the  velocity  of  small  particles 
that  are  assumed  to  move  with  the  flow.  If  these  particles  do  not  have  the  appropriate 
size  and  concentration  distribution,  the  experiment  results  may  be  adversely  affected. 
Particle  sizes  on  the  order  of  1  micron  are  often  quoted  as  adequate  for  most  gas  flows 
and  can  respond  to  turbulence  frequencies  exceeding  one  kHz.  Smaller  particles  may 


be  below  the  detectable  limit  of  the  instrument  but  will  contribute  to  the  background 
noise  while  larger  particles  that  are  easily  detectable  may  not  adequately  track  the  flow. 

It  w'as  determined  in  this  experimental  situation  that  seeding  of  the  flow  was 
not  absolutely  necessary.  Several  runs  were  made  with  no  seeding  supplied.  The  data 
obtained  in  these  runs  was  accurate  as  dust  particles  acted  as  the  flow-following  agent. 
However,  it  was  determined  that  by  using  a  flow  seeding  system  the  data  rate  was 
increased  by  a  factor  of  four  to  five.  The  system  used  consisted  of  a  paint  sprayer 
spraying  latex  particles  (on  the  order  of  3  microns)  suspended  in  alcohol. 

C.  COMPUTER-EXPERIMENT  INTEGRATION 

The  experiment  was  fully  integrated  with  a  Data  General  NOVA  minicomputer 
for  on-line  data  acquisition  and  analysis,  Figure  2.4.  A  main  program 
("SPOILER. FR")  was  written  to.  perform  all  tasks  and  record  all  data  required  during 
the  experiment  so  that  individual  tunnel  runs  proceeded  virtually  'hands  off.  The  main 
program  called  a  host  of  other  programs  during  each  run  in  order  to  accomplish 
various  tasks.  These  other  routines  are  listed  here: 

(1)  "SCANT. FR"  -  used  to  obtain  pressure  data  from  the  Scanivalvcs  and 
compute  free  stream  velocity  and  free  stream  Mach  number.  "SCANT. I  R" 
itself  calls  three  other  routines  in  order  to  determine  ambient  temperature, 
calibration  pressure  and  barometric  pressure: 

(a)  "TEMP.FR"  -  which  yields  ambient  temperature  from  a  thermocouple 
via  an  analog-to-digital  converter.  Ihe  thermocouple  is  located 
upstream  of  the  test  section. 

(b)  "CALPRHS.FR"  -  which  calculates  the  output  pressure  from  the 
calibration  vacuum  pump.  The  pressure  value  is  read  by  this 
program  using  an  analog-to-digital  converter.  I  his  value  is  then  used 
in  the  calibration  of  the 'Scanivalvcs. 

(e)  "NPARO.FR"  -  which  yields  the  value  for  current  atmospheric 
pressure  from  a  Paroscientilic  pressure  sensor  via  a  digital  I  O. 

(2)  "STD I’M. FR"  -  which  gives  the  current  time  and  date  in  order  that  the  data 
flies  being  created  by  the  main  program  can  be  labeled. 

(3)  "TRELAY.FR"  -  to  command  the  taskmaster  to  move  the  traverse  mechanism 
to  the  next  data  collection  location. 

(4)  "PLOTP.FR"  -  to  provide  on-line  data  analysis  by  allowing  any  of  five  curves 
to  be  plotted  immediately  after  a  given  run.  ' 

D.  EXPERIMENTAL  CONDITIONS 

T  he  free  stream  velocity  was  maintained  at  a  nominal  value  of  70  m  sec, 
corresponding  to  a  Reynolds  number  of  7.4  x  10"  based  on  the  airfoil  chord. 
Experiments  were  carried  out  with  the  spoiler  deflected  30  degrees  and  the  angle  of 
attack  at  0  and  .3  degrees.  Locations  where  data  was  obtained  are  shown  in  Figure  2.6. 


Figure  2.6  Data  collection  locations. 


In  order  to  determine  the  velocity  profiles,  data  was  acquired  at  varying  steps 
based  on  the  velocity  gradient.  Fine  step  sizes  were  taken  near  the  airfoil  for  boundary 
layer  profiles  and  near  free  shear  layers  in  the  wake  measurements.  Typically,  15  data 
points  were  taken  in  the  boundary  layer  and  20  points  in  the  wake.  The  4-beam  LDV 
arrangement  precluded  obtaining  data  closer  than  0.1  inch  from  the  airfoil  surface, 
however. 

The  boundary  layer  traverse  was  along  the  vertical  axis  of  the  test  section,  so  the 
traverse  direction  is  not  normal  to  the  airfoil  surface  as  are  usual  boundary  layer 
profiles.  Angle  of  attack  and  model  surface  curvature  make  the  measurement  axis 
incline  to  the  model. 


Ill.  DATA  ACQUISITION 


A.  FREE  STREAM  VELOCITY 

At  the  beginning  of  every  run,  the  main  experiment-driving  program 
("SPOILER.FR")  calls  for  the  Scanivalve  pressures,  ambient  temperature,  calibration 
pressure  and  atmospheric  pressure  to  be  read.  Assuming  isentropic  flow  in  the  test 
section  upstream  of  the  airfoil  and  away  from  any  shock  waves,  the  free  stream  Mach 
number,  M,  can  be  calculated  from  the  stagnation  pressure,  PQ,  and  the  static 
pressure,  P,  along  the  wall: 


M  =  {((P0/P)<^-1)/Y  -  1)(2/(y-I))}1/2 


For  y  =  1.4: 


M  =  {((P0,'P)(2/7)-  1)X5}1/2 


Then,  velocity  is  given  by: 


V  =  Ma 


where 


a  =  JyTt 


T  is  obtained  from  the  total  temperature,  T0,  using  the  following  formula: 


T  =  T  '{I  +  ((y-1 )  2)XM2} 
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Tq  is  measured  by  a  thermocouple  at  the  wind  tunnel  inlet. 

The  above  method  is  used  to  obtain  free  stream  Mach  number  and  free  stream 
velocity  for  every  run.  The  value  for  free  stream  velocity  is  used  in  the  data  reduction 
process  to  normalize  the  airfoil  velocities  obtained  from  the  LDV. 

B.  AIRFOIL  VELOCITIES 

Processing  of  the  Laser  Doppler  Velocimetcr  data  obtained  during  this 
experiment  was  carried  out  using  the  same  Data  General  Nova  digital  computer  and 
driving  program  used  for  all  other  experiment  functions.  "SPOILER. FR",  the  main 
program  for  this  experiment,  calls  a  digital  l  O  driver  program  ("DI024")  to  read 
words  as  fast  as  possible  from  two  cards  (these  cards  receive  the  digital  output  from 
the  two  Macrodyne  channels,  one  for  each  LDV  component).  As  soon  as  1000 
samples  (per  component)  have  been  taken  by  the  Macrodyne  processor  and  read  by 
DI024,  the  binary  output  from  each  sample  is  converted  to  decimal  using  a  decoding 
subroutine  ("DECODE l.FR"). 

As  mentioned  earlier,  the  Macrodyne  signal  processor  used  in  this  experiment 
utilizes  a  1  GHz  counter.  The  number  of  "counts"  required  for  the  passage  of  8  cycles 
is  sent  to  the  digital  I/O  in  the  form  of  a  10-bit  mantissa  and  4-bit  exponent,  namely 

Dg,  D8,...D0  X  2exP 

For  example,  if  400  "counts"  were  required  for  a  particular  Doppler  burst,  the 
output  would  be 


0000110010  0011 
Mantissa  EXP 

which,  upon  decoding,  reduces  to 

50  x  23  =  400 
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Since  the  clock  counts  are  produced  at  4  GHz  (the  1  GHz  clock  is  divided  into  4 
phases),  each  clock  count  is  0.25  nanoseconds.  Thus, 

PYn  (0.25) 

t_  =  Dg,  Dj,,...D(.  x  2exP  x  : - i  ns 

s  9  8  0  (8.0) 

is  the  period  of  the  signal  when  8  cycles  arc  measured. 

Converting  this  to  signal  frequency ,  fs,  yields 

1  25‘exP  x  1000 

f  =  —  =  -  MHz 

ts  Mantissa 

Since  the  signal  frequency,  fs,  is  composed  of  the  Doppler  frequency,  the  shift 
frequency  and  the  mixing  frequency,  the  Doppler  frequency  can  be  obtained  from 

^D  ~  *s  ^shift  +  *mix 

Once  the  Doppler  frequency  of  a  particle  is  known,  its  velocity  can  be  determined 

from 

velocity  =  const  x  f^ 

where  "const"  is  determined  based  on  the  fringe  spacing  or  direct  calibration. 

With  the  1000  samples  from  each  component  converted  to  decimal,  the  data  can 
be  reduced  to  average  and  rms  velocities  using  the  formulas  found  in  Table  II.  This 
data  (both  raw  and  reduced)  is  then  stored  for  future  analysis. 

While  most  data  analysis  is  performed  after  completion  of  the  experiment,  several 
on-line  analysis  routines  are  available  for  use  as  real  time  performance  checks.  For 
example,  1IISTOGRAM.FR"  is  a  program  available  to  first  read  a  given  number  of 
samples  from  one  of  the  Macrodyne  channels  and  plot  frequency  versus  the  number  of 


samples  at  each  frequency  in  histogram  fashion.  With  this  program  it  is  easy  to 
determine  whether  or  not  the  data  is  conforming  to  the  Normal  Distribution  (to  which 
ideal  data  should  conform).  Another  on-line  analysis  tool  which  is  built  into  the  main 
program  "SPOILER. FR"  is  the  subprogram  "PLOTP.FR".  This  routine  gives  the 
experimenter  the  option  to  have  any  of  five  velocity  profiles  (i.c.,  L'ave  Uqq, 
Vave/Voo>  etc>)  plotted  on  the  screen  immediately  after  completing  an  experimental 
run.  This  subprogram  was  used  extensively  to  immediately  inspect  the  quality  of  each 


TABLE  II 

VELOCITY  PROFILE  AND  TURBULENT  STRESS  EQUATIONS 
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C.  REYNOLDS  SHEAR  STRESSES  AND  TURBULENT  STRESSES 

Table  II  shows  the  equations  used  to  determine  mean  velocities,  turbulent 
stresses  and  Reynolds  shear  stress  throughout  this  experiment.  The  value  for  N  was 
generally  1000  since  this  was  the  number  of  velocity  samples  taken  at  each 
measurement  location  (for  each  component). 


IV.  RESULTS  AND  DISCUSSION 

A.  INTRODUCTION 

One  of  the  purposes  of  this  experiment  was  to  demonstrate  the  benefits  of  the 
two-component  Laser  Doppler  Velocimeter  by  employing  this  type  of  system  to 
measure  the  flow  field  of  an  airfoil  with  a  deflected  spoiler.  In  order  to  present  the 
results  of  these  measurements  effectively,  a  graphics  program  called  "PLOT-3D"  was 
used  to  plot  the  data  obtained.  "PLOT-3D"  is  available  on  the  VAX  system  at  the 
NASA  Ames  Research  Center. 

A  data  interpolation  routine  written  by  Chvang  Shcng  Lee  of  Stanford  University 
was  also  used  to  interpolate  the  measured  data  into  a  grid  system  so  that  it  could  be 
plotted.  All  of  the  contour  plots  presented  use  the  interpolation  routine  to  interpolate 
in  both  the  x  and  y  directions.  However,  some  of  the  velocity  vector  plots  use 
interpolation  in  the  v  direction  only.  This  was  done  because,  in  some  cases,  linear 
interpolation  of  the  data  gives  a  false  impression  of  what  the  How  field  is  actually 
doing. 

In  all  cases,  the  data  presented  in  this  section  is  for  a  =  0  0  and  6  =  30°. 

B.  FLOW  FIELD  ANALYSIS 

Velocity  measurements  were  taken  at  a  variety  of  points  along  the  upper  surface 
ahead  of  and  behind  the  spoiler.  Figure  4.1  shows  a  plot  of  the  data  obtained  over  the 
rear  portion  of  the  airfoil.  Note  the  low- velocity,  recirculating  llow  region  directly 
behind  the  spoiler. 

While  the  data  for  Figure  4.1  was  interpolated  in  the  v  direction  only,  the  data 
for  Figure  4.2  was  interpolated  in  both  the  x  and  y  directions.  It  can  be  seen  in  Figure 
4.2  that  interpolation  in  the  x  direction  yields  some  fairly  unrealistic  velocity  vectors 
due  to  the  relatively  low  number  of  measurement  locations  as  compared  to  the  number 
of  grid  points.  For  example,  those  directly  ahead  of  the  spoiler.  On  the  other  hand, 
complete  interpolation  (in  both  the  x  and  y  directions)  docs  serve  to  paint  a  good 
overall  picture  of  what  the  llow  field  is  doing  based  on  a  finite  number  of  velocity 
measurements. 
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C.  WAKE 


The  wake  behind  the  spoiler  was  measured  with  the  two-component  LDV  system 
at  numerous  locations  to  determine  the  mean  and  rms  profiles.  The  results  from  the 
wake  are  presented  in  this  section. 

1.  Wake  Profiles 

The  mean  longitudinal  and  transverse  wake  profiles  measured  by  the  LDV  for 
a  =  0°  and  6  =  30°  are  presented  in  Figures  4.3  and  4.4.  The  symbols  are  the 
measured  data  and  the  lines  represent  the  best  fit  curves  of  the  data.  As  expected,  the 
width  of  the  wake,  the  velocity  defect  and  transverse  velocity  all  decrease  in  the 
downstream  direction. 

The  Reynolds  shear  stress,  longitudinal  and  transverse  turbulent  stresses 
measured  by  the  LDV  are  presented  in  Figures  4.5  to  4.7.  Obvious  in  all  three  of  these 
figures  is  the  fact  that  the  fluctuating  region  appears  to  grow  from  the  x,  C=  1.0 
(trailing  edge)  location  to  x/C=  1.25  then  decrease  as  it  approaches  x/C=  1.5.  The  "S" 
shape  profile  in  the  Reynolds  shear  stress  is  typical  for  a  two-dimensional  wake 
[Ref.  l:p.  29J.  The  longitudinal  and  transverse  turbulent  stresses  show  two  peaks  at  the 
locations  where  the  two  free  shear  layers  (corresponding  to  the  airfoil  and  the  spoiler) 
oscillate  in  the  near  wake. 

Figures  4.8  and  4.9  show  the  velocity  vectors  in  the  wake.  A  scaled  profile  of 
the  airfoil  is  superimposed  in  each  figure  for  reference.  The  data  in  Figure  4.8  was 
interpolated  in  the  y  direction  only.  Figure  4.9  was  interpolated  in  both  directions  and 
depicts  the  mean  wake  velocity  profile  quite  well. 

2.  Contour  Plots 

f  igures  4.10  to  4.16  depict  "contour"  plots  of  various  quantities  in  the  airfoil 
wake.  Figure  4.10  shows  regions  of  equivalent  L'-magnitude  (U  is  referenced  to  free 
stream  velocity).  Note  that  there  are  regions  where  the  ratio  U  U  x>  's  greater  than 
unity.  T  his  is  due  to  tunnel  blockage  effects  which  cause  the  How  to  accelerate  more 
than  normal  around  the  airfoil. 

Figure  4.11  depicts  the  regions  of  equivalent  V-magnitudc  (also  referenced  to 
free  stream  velocity).  A  minus  sign  on  a  contour  line  in  Figure  4.1 1  implies  that  the  V 
component  of  velocity  is  negative  (downward)  in  that  region. 

In  Figure  4.12,  the  total  velocity  magnitudes  in  the  wake  are  contoured. 
Again,  contour  lines  with  values  greater  than  one  arise  due  to  blockage  effects. 


x/C 


Figure  4.9  Velocity  vectors  in  the  wake 
(interpolation  in'  both  directions). 


Figure  4.15  Contour  plot  of  transverse  turbulent  stresses 

in  the  wake. 


Figure  4.13  shows  how  rapidly  the  flow  angles  change  in  the  near  wake  (the 
closer  the  contour  lines,  the  steeper  the  gradient). 

Longitudinal  and  transverse  turbulence  and  Reynolds  shear  stress  are 
contoured  in  Figures  4.14  to  4.16.  Figure  4.14  clearly  shows  the  two  regions  where  the 
turbulence  intensity  is  the  highest.  These  two  regions  are  probably  caused  by 
counter-rotating  vortices  originating  at  the  spoiler  trailing  edge  (clockwise  vortex)  and 
the  airfoil  trailing  edge  (counterclockwise  vortex). 

D.  COMPARISON  WITH  HOT-WIRE  ANEMOMETER  DATA 

One  purpose  of  this  thesis  was  to  compare  the  data  obtained  in  this  experiment 
with  the  hot-wire  data  obtained  under  similar  conditions  by  Lee  [Ref.  5].  Figures  4.17 
to  4.21  are  LDV/hot-wire  comparisons  of  longitudinal  and  transverse  velocity, 
longitudinal  and  transverse  turbulent  stress,  and  Reynolds  shear  stress  at  three 
locations  in  the  wake. 

It  is  important  to  mention  that,  in  the  previous  experiment,  no  hot-wire  data  was 
obtained  in  the  region  of  reverse  flow  between  x/C=1.0  and  x/C=1.25.  Hot-wire 
anemometry  is  virtually  unusable  in  this  region  due  to  the  fact  that  hot-wires  rectify  the 
signals,  which  include  the  instantaneous  forward  velocity  as  well  as  the  reverse  velocity. 
Consequently,  velocities  measured  by  hot-wires  in  reverse  flow  regions  appear  higher 
than  they  actually  are. 

Except  for  minor  grid  alignment  errors  in  the  y  direction,  most  of  the  curves  are 
in  very  good  agreement.  As  expected  however,  the  longitudinal  turbulence  stress 
curves  at  x  C  =  1.25  (Figure  4.19)  and  the  transverse  turbulence  stress  curves  at 
x  C=  1.25  (Figure  4.20)  show  considerable  disagreement.  In  the  longitudinal 
turbulence  case,  the  hot-wire  data  shows  more  turbulence  in  the  middle  of  the  wake 
while  in  the  transverse  turbulence  case,  the  hot-wire  data  shows  less  turbulence  in  the 
middle  of  the  wake.  This  is  again  due  to  the  fact  that  hot-wire  measurements  give  a 
misleading  velocity  profile  close  to  the  trailing  edge  because  of  the  reverse  How  in  that 
region  [Ref.  l:p.  70). 

The  hot-wire  anemometer  is  quite  accurate  in  regions  where  there  is  no  reverse 
flow.  However,  in  order  to  obtain  reliable  data  on  the  mean  velocities  and  turbulent 
stresses  in  the  reverse  (low  regions  and  the  near  wake  of  this  airfoil-spoiler  system. 
Laser  Doppler  Yelocimetrv  yields  better  results. 


0.12 


2  0 


.25  x/C=1.5  x/C=2.0 


.'am.v-*av"- ""vV-TT-^TW 's -.  -T-.'q'vv.vvA 


E.  ERROR  ANALYSIS 

In  this  experiment,  1000  velocity  samples  were  taken  by  the  LDV  in  each 
direction  (U  and  V)  at  every  measurement  location.  The  purpose  in  taking  so  many 
samples  was  to  provide  enough  values  in  the  fluctuating  regions  to  ensure  that  an 
accurate  statistical  mean  could  be  computed.  Consequently,  in  the  data  reduction 
process  all  1000  values  were  averaged  to  obtain  values  for  L'ave  and  Vave  at  each 
location.  The  equations  in  Table  II  were  then  used  to  obtain  Reynolds  shear  stress, 
longitudinal  and  transverse  turbulent  stress  at  these  locations  (still  using  all  1000 
samples). 

To  determine  its  effect,  a  different  data  processing  method  was  tested  on  the  raw- 
data  from  several  runs.  In  this  procedure,  only  the  samples  which  fell  within  one 
standard  deviation  of  the  mean  were  considered  (approximately  67%  of  the  1000 
samples)  in  computing  L'ave,  Vaye,  L'rms  and  Vrms.  Again  the  equations  in  Table  II 
were  employed,  just  as  before.  This  new  process  had  little  effect  on  L'e  but  Vave 
changed  considerably  (from  -1.3  to  -11.0  at  one  location  tested).  However,  the  V 
values  did  not  always  decrease  (one  location  changed  from  8.7  to  16.9).  Therefore,  the 
changes  were  not  due  to  the  removal  of  a  bias  error.  The  longitudinal  and  transverse 
turbulence  stress  values  decreased  in  all  cases.  This  was  expected  since  the  biggest 
contributors  to  these  values  are  the  velocity  terms  which  dilfer  greatly  from  the  mean. 

This  new  data  processing  procedure  did  not  improve  the  data  quality  by 
removing  any  biases  (based  on  the  fact  that  some  values  of  L'avc  and  Vavc  increased 
while  others  decreased  on  using  the  new  process).  As  a  matter  of  fact,  this  process  had 
a  negative  effect  on  the  data  quality  in  that  it  neglects  velocity  extremes  which,  in  turn, 
lowers  the  apparent  turbulence  intensity  of  the  flow  field.  For  this  reason,  the  original 
data  processing  procedure  (consider  all  1000  samples)  was  used  in  this  experiment. 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

A  two-component  LDV  system  was  set  up  for  the  purpose  of  measuring  mean 
and  fluctuating  quantities  of  an  airfoil  with  a  deflected  spoiler.  The  results  of  this 
experiment  have  shown  that  this  system  provides  an  efficient  and  effective  means  of 
flow  field  measurement,  including  regions  of  reverse  flow  and  large  velocity 
fluctuations.  Comparison  with  hot-wire  data  demonstrates  the  LDV's  superiority  over 
the  hot-wire  in  measuring  in  the  near  wake  region. 

The  results  also  clearly  show  two  recirculating  regions  in  the  wake  of  the  airfoil. 
One  vortex  originates  at  the  trailing  edge  of  the  airfoil  itself  while  the  other  originates 
at  the  trailing  edge  of  the  spoiler.  The  turbulence  region  is  concentrated  in  the  near 
wake  and  dissipates  rapidly  in  the  downstream  direction. 

B.  RECOMMENDATIONS 

In  order  to  establish  a  more  complete  data  base,  it  is  recommended  that  more 
LDV  measurements  be  obtained  at  other  than  a  =  0°  and  6  =  30°.  These  results  are 
important  to  the  modeling  of  this  flow  field.  Also,  to  remove  the  possibility  of  tunnel 
blockage  errors,  it  is  necessary  to  perform  further  studies  in  a  large  scale  wind  tunnel. 

When  using  an  LDV  to  measure  a  turbulent  flow  field  such  as  in  this  experiment, 
a  large  number  of  samples  should  be  taken  at  each  location.  Furthermore,  unless 
attempting  to  remove  a  velocity  bias  all  of  the  samples  taken  should  be  considered  in 
computing  L'ave,  Vave,  etc. 
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APPENDIX 

COMPUTER  PROGRAMS 


PURPOSE  -  (1)  TO  CONVERT  A  24  HOUR  CLOCK  TO  A  12  HOUR  CLOCK 
-  (2)  TO  ADD  A.M.  OH  P.tl.  TO  THE  “ADJUSTED"  TIME. 


SUBROUTINE  STDTM  ( I  DATE, ITIME, IVAL) 
DIMENSION  I  DATE ( 3 ) 

DIMENSION  IT  1 HE!  3  ) 

CALL  DATE  ( 1  DATE . I  ERROR ) 

CALL  TIME  ( ITIME, I  ERROR) 

IF  ( IT1ME( 1 ) -GE. 13)  CO  TO  30 
IF  ( ITIME! 1 ) .10. 12)  CO  10  10 
CO  TO  20 

10  IF  ( ITIMEI2) .GT.OO)  CO  TO  30 
IF  ( 1TIMEI3) .Cr.OO)  CO  TO  30 
20  IVAL  =  "AM” 

CO  TO  40 
30  IVAL  =  "PM" 

IF  ( ITIME! 1 ) .LT. 13)  CO  TO  40 
I T I  ME ( 1 )  =  ITIME! 1 )  -  12 
CO  TO  50 

40  IF  ( ITIME! 1 ) .CE. 01 )  CO  TO  50 
ITIME! 1)  =  ITIME! 1)  +  12 
50  CONTINUE 
RETURN 
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*****  ************************* 
************************:;••***** 


**  ** 

**  SPOU.ER.fr  ** 

**  ** 

**  MIKE  FOREMAN  ** 

** 


*************************<,:**:|;£ 

*  ***  S**  *****  ********  ***  ******* 

LASER  DATA  ACQUISITION  PROGRAM 

RLDll  SPOILER  STDTM  D1024  6 ELI  BO 

DIMENSION  I  FI LE( 4 ) ,  IDT(IO),  VAR(IO) 

DIMENSION  ITIMEtO),  1  RATH  3 ) 

DIMENSION  1  RAY  I ( 1 000 ) ,  IRAY2II000) 

DIMENSION  1 1)  I  N  I  (16)  ,  1RIN2U6),  VEL1  (1 000  )  .  VEL2  (  1  000 ) 

D 1  MENS  ION  V  (  0 . 20  )  ,  X  (  3</ )  ,  Y  <  30 ) 

INITIALIZE 

DO  10  1=1,3 
I  DTI  I ) =0 
VAR  I  I )  =  1 0 . 

10  CONTINUE 

ACCEPT  “RUN  NUTIRER  (ZERO  TO  QUIT):  “,IRUN 
IF<  IRUN.I.E.O)  CO  TO  009 

15  ACCEPT  " NUMBER  OF  SAMPLES  PER  COMPONENT:  *  , NTJM 
I  E  (  NIJM .  I.F. .  1 000 .  AND .  HUM .  GE .  0  )  GO  TO  20 
TYPE  "NUMBER  OF  SAMPLES  MUST  UE  BETWEEN  1  AND  1000“ 

CO  TO  15 

20  ACCEPT  "MIXING  FREQUENCY  *1:  “ .FREOI 
ACCEPT  "MIXING  FREQUENCY  *2:  “ , FRE02 

CONI =16.65 
C0N2= 16.65 
VAR ( 1 ) =CON 1 
VAR12 ) =C0H2 
C 

C  SURDEVICE  ADDRESSES 
0 

1SUBI =46K 
I SUB2=47K 
C 

C  TIME  PROCILAH 
C 

C  ALL  STDTM ( I  D  ATE , IT IKE , I VAL ) 

C 

DO  30  1=1,3 
I  DTI  I )  =  I  DATE I  I ) 

I  DTI  I+3)=HTIiE(l  ) 

30  CONTINUE 

I  DTI  7  )  =  I VAI. 

I  DTI 3)  =  I  RUN 
V  All  I  0  )  =  FREQ  1 
VAR I  1 0 ) =  FREQ2 
C 

35  TYPE  "WILE  THIS  DE  A  :  II)  WAKE  PROFILE" 

ACCEPT  *  12)  BOUNDARY  LAYER  PROFILE  • . I  PROF 

G 

I F I  1  PROF . EQ . 1 . OR . 1  PROF . EO . 2 )  GO  TO  40 
TYPE  "MUST  ENTER  1  OR  2" 
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CO  TO  35 

40  IK< 1PROF.EQ.2)  CO  TO  41 

ACCEPT  "X-LOCATI ON  FOR  THIS  RUN  (XT/C):  -.VARO) 

CO  TO  42 

41  ACCEPT  "X-LOCATI ON  FOR  THIS  RIJN  (X'C):  ",VAR(3) 

42  ACCEPT  "NUMBER  OF  Y-LOCAT 1  OHS :  ",NMAX 
1F( I  PROF. EO. 2)  CO  TO  45 

ACCEPT “STANDARD  WAKE  PROFILE  Y-LOCATIONS?  *1=YES,0=NO)  " , LL 
I F  <  LI, .  F.Q .  1  )  CO  TO  44 

TYPE" I NPUT  Y-LOCATIONS  IN  INCHES  ( LOWEST  POINT  FIRST):" 

DO  43  1  =  1,  NfLAX 
WRITE  (10,49)  I 
ACCEPT"  “ , Y< I ) 

43  CONTINUE 
CO  TO  50 

44  Y ( 1 )=-l .0 
Y ( 2 ) = -0 . 8 

Y  *  3 ) = -0 . 6 
Y ( 4 ) = -0 . 46 

Y  <  5 ) = -0 . 34 

Y  ( 6 )  =  -0 . 2 
Y*7) =-0 . 1 
Y<  8 ) =0 . 0 
Y<  9 ) =0 . I 
YC 10)=  .2 
Y* 1 1 )= .3 

Y  < 12)= .4 
Y< 13)= .5 
Y ( 14)= .6 
Y< 15)= .7 
Y< 16)= .0 
Y(  17)  =  .9 
Y< 1 8 ) = 1 .0 
Y ( 1 9 )  =  1  .2 
Y ( 20  )  =  1  .4 
CO  TO  50 

45  TYPE  STANDARD  BOUNDARY  LAYER  PROFILE  Y-LOCATIONS?  " 

ACCEPT"  < 1 =  YES ,  0  =  NO>:  ",LK 

IFO.K.EQ.  1  )C0  TO  47 

TYPE" INPUT  Y-LOCATIONS  IN  INCHES  (LOWEST  POINT  FIRST):" 

DO  46  1  =  1  ,  NflAX 
WRITE* 10,49)  I 
ACCEPI . ,Y(  I  ) 

46  CONTINUE 
CO  TO  50 

47  Y< 1 ) = . 1 
Y(2)  =  .2 
Y(3)=  .3 
Y<4)=  .4 
Y<  5  >  =  .5 
Y ( 6 ) =  .6 
Y ( 7  )  =  .7 
Y<8) = .8 
Y<9)= .9 

Y  (  1 0  )  =  1  . 0 

Y  (  1  I  )  =  1  .  1 
Y( 12) = I .2 
Y ( I  3 )  =  1  .3 
Y( 14)= I .4 
Y< 1 5 ) = 1 .5 

49  FORM \T(6X,"Y(",I2,"):“) 

50  W|UTE(  10, 125) 

125  FORMAT*"  "."BEGIN  (SCANI.FR)  >>",Z) 

CALL  FSWAP * "SCAN! .SV" ) 
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ACCEPT  “DO  YOU  WANT  RAW  DATA  STORED?  <1=YES,  0  =  NO)  “,IRAW 
I DT ( 9 ) = I RAW 
I DT ( 1 0 )  =  I  PROF 

BECIN  TAKINC  DATA 


N  =  0 

57  CALL  DIOI I ( ISUDI , IRAY1 , IRAY2.NUM, 1ER) 
DECODE 
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65 


70 


75 


AVE1 =0. 

AVE2=0 . 

DO  90  1=1. NUM 
DO  65  J=I , 16 

1B1N1 (J  >  =  I AND< IRAY1 ( I > ,2**< J-l  ) ) 
I  B I  N2 <  J  )  =  I ANIH I RAY2 ( I ) ,2** (J-l ) ) 
IFUBINI  (J)  .NE.O)  1  D I  N I  (.I>  =  1 
1KIBIN2CJ)  .NE.O)  1B1N2(J)=1 
CONTINUE 
I  BN  1 =0 
1 BN2=0 

DO  70  J=0,9 

IBNI  =  1  BN  I  + 1 B1 N 1 ( J+ 1 )*2**J 
1BN2=IBN2+1BIN2(J+1 )*2**J 
CONTINUE 
1  EXP  I =0 
I EXP2=0 

DO  75  J=10, 13 

I  EXP 1=IEXP1+1BIN1 (J  +  l  )*2**<J-10) 
I EXP2= 1 LXP2  + I B1 N2 ( J+ 1 )*2**<J-10> 
CONTINUE 


VEL1 ( I )=C0N1*(2.**(5-IEXP1 ) *  1 000 . /FLOAT ( IBNI J-40.+FREOI ) 
VEL2 ( I )=CON2*<2.**(5-IEXP2)*1O0O./FLOAT< lBN2)-40.+FREQ2) 
AVE1 =AVEI +VEL1 ( I ) 

AVE2= AVE2+VEL2( I ) 

90  CONTINUE 

AVE 1 =AVE1/NUM 
AVE2= AVE2/NUM 

CHECK  VELOCITIES  ARE  IN  LIMITS 


NN  =  N+ 1 

WRITE( 10,92) NN , AVE2 , AVE 1 

92  F0RMAT< IX. “LOCATION:  " . 1 2 , 4X, *U( AVE) =  * , F9 . 5 ,4X, • V ( AVE) =  ■ 
I F ( AVE 1 . CT . 35 . )  GO  TO  95 
IF(AVE2.CT. 140. )  CO  TO  95 
CO  TO  90 

95  WHITE  FHEE(IO)  "BAD  DATA:  LOCATION  " ,NN 

ACCEPT “RETRY  OR  SKIP  TO  NEXT  LOCATION?  ( 1 =RETRY,  0=SKIP>: 

I F ( I  TRY . EQ . 0 ) CO  TO  90 
CALL  DIORM 
CO  TO  57 

90  I F < IRAW.NE. 1 )  GO  TO  100 

STORE  RAW  DATA 

CALL  FOPF.N (5.  "COUNTER"  ) 

READ  BINARY<3)  NJ 


,F9.5) 


" , I TRY 


CALL  FCLOSE( 5 ) 

NJ=NJ+1 

CALL  FOPEN <5, "COUNTER" ) 

WRITE  BINARY ( 5 )  NJ 
CALL  FCLOSE( 5 ) 

NFI LE-NJ 

IF(N.EO.O)  JFILE=NFILE 
I  R=MOD( NF I LE ,  10) 

1 S=NFI LE/1O0 
I L  = (NFI LE- I S*  1 00 )/l 0 
I  F 1  LE (  1 )=’RA’ 

I F 1 LE ( 2 ) = ' WD' 

IFILE<3)= ( IS+601O*4OOK+( IL+60K) 

I F 1 LE  1 4  >  = ( 1R+60K)*400K 
C 

CALL  FOPEN ( 5 , I  FI LE) 

WRITE  BINARY(5)  ( I DT( I ) , I = I . 10 > 

WHITE  BI NAHY ( 5 )  ( VAH( 11,1=1,10) , NUM , Y(NN ) 

WHITE  BINARY (5)  ( 1  RAY  1 ( I ) . 1  =  1 . NUM ) 

WHITE  BI  NARY (  5  )  (  1 1LAY2 (  I  )  ,  I  =  1  ,  NUM ) 

CALL  FCLOSEt 5 ) 

C 

100  WHITE  (10,99) 

99  F()IlfL\TC  ","***  BEGIN  RMS  CALCULATIONS  >>*,Z) 

urms=o. 

vitns  -  o . 

UVRMS=0. 

DO  110  1=1, NUM 

URMS=URMS+(  ( VEL2  (  I  )  )  **2- ( AVE2  )  **2  ) 

VRMS= VRMS+ (  ( V  EL  I  (  1  )  **2-  ( AVE1  >**2)  ) 
UVIU;S=UVRM3+<VELl ( 1 >*VEL2( 1 ) ) - ( AVE1 *AVE2 ) 

110  CONTINUE 

URi;S=  (URMS/FLOAT(NUN)  ) 

VHMS  =  ( V1U1S/FL0AT(NUM)  ) 

UV  It  I  IS  =  ( UVRMS/ FLOAT  (  NUM )  ) 

WHITE  (10,111) 

111  F()ilM\T(  "  +  END  RMS  CALCULATIONS") 

N  =  N  +  1 

V( 1  ,  N  )  =  AVE2 
V ( 2 , N ) =  AVE1 
V(3,N)=URMS 
V  ( 4  ,  N  )  =  VIUIS 
V ( 3 , N ) =U VRMS 
IF(N.CE.NHAX)  GO  TO  120 
C 

C  MOVE  TRAVERSE 
C 

CALL  FSWAP ( "TRELAY. SV " ) 

ACCEPT "TRAVERSE  MOVED... ANY  NUMBER  TO  CONTINUE.  * 

0 

CALL  1)1  ORIl 
CO  TO  57 
C 

C  OBTA I N  TEMP ,  STAG  PRES . .  MACH  * ,  ETC 
G 

120  I F ( IRUN.EQ. I COO)  GOTO  127 
C A LL  FSWA P ( " TRELAY . SV "  ) 

WRITE ( 10,126) 

126  FORMAT! “TRAVERSE  RESET  TO  INITIAL  POSITION") 

127  CALL  FOPEN (5, "SCAN  I  "  ) 


READ  It  1  NAHY  <  5  )  VAR  ( 4  )  ,  VAR  (  5  )  ,  VAR  ( 6 )  ,  VAR  ( 7  )  .  VAR(  8 ) 


CUE  rCI.0SE<5) 


C 


G  MUTE  DATA  INK)  PERMANENT  DATA  FILE 
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CALL  F0PEN(5,  'TERMCOUNT"  ) 

READ  Bl NARY ( 5 )  ME 
CALL  FCLOSE(5) 

MUTE  FREE(IO)  "LAST  DATA  FILE  USED:  *  ,MF 
MUTE  FREE(IO)  "HUN  NUMBER :  ".IRUN 

ACCEPT  “'DATA'  FILE  NUMBER  (NEC  *  TO  SKIP):  " ,MFILE 
IF(MFILE.LT.O)  CO  TO  599 
CALL  FOPEN  (  5  ,  "PEIUICOUNT"  ) 

MUTE  B1NARY(5)  Ml  lLE 
CALL  FCLOSE<  5 ) 

IR=MOD(MFILE,  10) 

IL-riFlI.E/10 
1 F I LE  < 1 )  =  ’  DA ’ 

I F 1 LE( 2 ) = ’ TA* 

I  F I  LE ( 3 )  =  (  I  L+60K ) *4Q0K+  (  IR+60K) 

I F I LE ( 4 )  =  0 

CALL  FOPEN (5. I  FILE) 

MUTE  BINARY!  5)  (  I DT (  I  )  ,  I  =  I  ,  10 ) 

MUTE  Bl  NARY ( 5 )  ( VA1K  I  )  .  I  =  I  .  1 0 )  . NUM. NMAX .  JFI LE,  NFI  LE 
MUTE  BINARY ( 5 )  <  ( V ( I  ,  J  )  .  J  =  I  . NMAX)  ,  I  =  1  . 5 ) 

WRITE  BINARY! 5)  ( Y ( I ) , I = I , NMAX ) 

CALL  EC LOSE! 5 ) 

599  WRITE  FREE!  10)  “AMBIENT  TEMP:  “  ,VAR(4)  ,  "  DECF" 

MUTE  FREE!  10)  "ATMOSPIIER  1C  PRESSURE:  “  ,VAR(5)  ,  "  PSI" 

MUTE  FREE!  10)  "Sl'ACNATIOU  PRESSURE:  ",VAR(6),"  PSI" 

WRITE  FREE! 10)  “FREE  STREAM  MACH  NUMBER:  " ,VAR<7) 

MUTE  FREE  (10)  "FREE  STREAM  VELOCITY:  “ .VAR(8),"  M'S” 

DO  630  1=1, NMAX 
Y ( 1  )=Y( I )/6. 

650  CONTINUE 

ACCEPT "DO  YOU  DESIRE  A  CURVE?  (1=YES,  0=NO) :  ".ICURVE 
IF! 1  CURVE. EO.O)  CO  TO  973 
MUTE  FREE!  10) 

600  TYPE" WHAT  PROFILE  DO  YOU  DESIRE?" 

TYPE"  (1)  U(AVE)/U( INFINITY)" 

TYPE"  (2)  V(AVE)/U( INFINITY)" 

TYPE"  (3)  U<RHS)/<U(INFINITY>**2)" 

TYPE"  (4)  V(RMS)/(U! I NFI NITY)**2) " 

ACCEPT"  (3)  UV(1UIS)/(U( I NF I N ITY > **2 )  " . ISEL 

BIC1U0. 

DO  7uO  1=1, NMAX 

1  F(  1SEL.LE.2)  X!  I  )=V<  ISEL,  I  >/VA!UB> 

IF! ISEL.CT.2)  X! I )=V( ISLE. I )/ ( VAR! 0)**2 ) 

IF(X( I ) .CT.BICX)  BICX=X(l) 

700  CONTINUE 

CALL  FOPEN (5, “PROFILE" ) 

WRITE  BINARY! 5)  NMAX , I  PROF , I SEL , Bl CX 
WRITE  BINARY ( 5 )  (X! I  )  ,  I  =  I  .NMAX) 

WRITE  BINARY! 5)  ( Y ( I ) , I = 1 , NMAX ) 

CALL  I  CLOSE! 5) 

CALL  I  SWAP ( “ PLOTP . SV “ ) 

ACCEPT'ANY  *  TO  CONTINUE  “ .ANY 
TYPE  "DO  YOU  DESIRE  ANOTHER  GRAPH?  " 

ACCEPT'  ( I =  YES ,  U  =  NO) :  ",ICR 

IF! I0R.EQ.1 )C0  TO  600 
973  WRITE  FREE! 10) "ENSURE  PRINTER  IS  ON" 

WR I TE  < 1 2 . 9G0 ) I  DATE , 1 T I  ME . I  V AL 

9U0  FORMAT! IOX, "DATE:  " , 2 ( 12,  'V  “  )  ,  1 4 , 5X, "TIME:  " , 2 ( 12, " : " ) , 1 2 , 1 X , A2 ) 
WRITE  FREE! 12) 

WRITE  FREE! 12) "RUN  NUMBER:  ",!RUN 
WRITE  FREE! 12) 

IF!  II'ROF.EO.  1  )  WRITE  FREE!  I  2 )  "WAKE  PROFILE" 
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IF( irnOF.LQ.2)  WRITE  FREE ! I  2 ) “ HOUND, MIY  LAYER  PROFILE" 

WRITE  FREE< 12) “  X-LOOATION:  ",VAH!3) 

WRITE  FREE( 12) 

WRITE  FREE! 12) “NUMBER  OF  Y-LOCATIONS  TAKEN:  " , NMAX 
WRITE  FREE! 12) 

WRITE! 12,901 ) 

981  F0RMAT12X,  “LOCATION"  ,  10X,  "  U  ( AVE )  “  .  I  OX ,  "V(AVE)  “  ,  10X,  "UIRMS) "  , 

81 IX, “VIRUS) "  ,  1  IX,  "UVIRNS)  “  ) 

DO  935  1=1. NMAX 

WRITE! 1 2 , 984 ) Y ( I),V!l,I),V!2.I).Vi:i,I),V(4,I),V!5,I) 

984  FORMAT < 2X , F8 . 6 , 9X . F9 . 4 , 8X , F9 . 4 . 7X , F 1 0 . 4 , 7X , F 1 0 . 4 . 7X . F 1 0 . 4 ) 

905  CONTINUE 

WRITE  FREE! 12) 

WRITE  FREE ! 1 2 ) *  NUMBER  OF  SAMPLES  PER  Y-LOCATION:  ",NUM 
WRITE  FREE!  12)  "fl!X  FREQ  •»  1  =  "  . FRF.Q  1  ,  "  NIX  FREQ  #2=  "  , FREQ2 
WRITE  FREE! 12) " VEL  CONST  *1=  ".CONI,"  VEL  CONST  •I-  " ,C0N2 
WRITE  FREE! 12) 

IF! IRAW.EQ.0)  WRITE  FREE! 12) “RAW  DATA  HAS  NOT  BEEN  STORED" 

IF! IHAW.EQ. I )  WRITE  FREE! 12) “RAW  DATA  IS  STORED  IN  'RAWD’  FILES 
SJFILE, *  THRU  “ .NFILE 

WRITE  FREE! 12) “PERMANENT  DATA  FILE:  DATA" ,MF ILE 
WRITE  FREE! 12) 

WRITE  FREE! 12) "AMBIENT  TEMP=  “,VAR!4),"  DECF " 

WRITE  FREE!  12)  "ATMOSPHERIC  PRF.SSURE=  M,VAR!5)."  l’Sl  " 

MUTE  FREE! 12) "STAGNATION  PRESSURE=  “,VAR16)("  PSI " 

WRITE  FREE!  12)  "FREE  STREW  M  MACII  *  =  *,VAR(7) 

WRITE  FREE! 12) "FREE  STREAM  VELOCITY=  H,VAR!0),"  M/S" 

999  STOP 
END 


63 


c  *************************** 

G  *************************** 

C  **  ** 

C  **  "SCANI.FR."  ** 

C  **  ** 

C  **  MIKE  FOREMAN  ** 

C  **  ** 

G  *************************** 

C  *************************** 

c 

C  "SCANI.FR"  TAKES  THE  DATA  FROM  THE  6 
C  SCAN  I  VALVES  AND  WRITES  TUIS  DATA  (IN  COUNTS) 

C  ON  THE  SCREEN  AND  INTO  A  DATA  FILE. 

C 

DIMENSION  IRAY<1200>,  Y<6,0:23> 

C 

C  TIIE  SCAN  I  VALVES  ARE  CONNECTED  TO  CHANNELS  1  TRROUCH  6: 

G 

N START  =  1 
LIMIT  =  6 
NUN  =  L I M I T-NSTART+ 1 
C 

C  ZEROIZE  THE  MATRIX 
C 

DO  10  1  =  1.6 
DO  10  J  =  0.20 
Y< I  .  J)  =  0 
10  CONTINUE 
C 

0  SCAN  I  VALVE  RELAY  INFORMATION 
C  I  DEV  ICE  =  DEVICE  NUMBER  (25K  FOR  SCAN  I  VALVE) 

C  RELAY  NUMBERS:  HONE  =  RELAY  *0 
C  STEP  =  RELAY  *2 

C  TO  CONVERT  RELAY  NUMBER  TO  DINARY  WORD  USE: 

C  N=2**(16-N> 

C 

I  DEV  ICE  =  23K 

N 1 =2  ; STEP  =  RELAY *2 

I  STEP  =  2**< 16-N1 ) 

N2  =  3  iHO.lE  =  RELAY  >-0 

I  HONE  =  2** ( 1 6-N2 ) 

C 

C  RESET  ALL  BITS  TO  ZERO 
G 

CALL  DOUTW ( I  DEV  I CE , 0 ) 

C 

C  SETUP  THE  A/D  DEVICE  WITH  THE  STARTINC  AND  ENDING  CHANNELS 
C 

CALL  FADSKNSTART,  LIMIT,  I  Ell) 

C 

C  RESET  THE  SCAN  I  VALVE  TO  THE  HONE  POSITION 
C  (ENERCIZE  THE  RESPECTIVE  RELAY,  THEN  DEENERCIZE  IT) 

C 

CALL  DOUTW( IDEVICE, inONE) 

C  ALL  FDELYI  5) 

CALL  DOUTW ( I  DEVICE,©) 

G 

C  TAKE  200  SAMPLES  AT  EACH  PORT 


DO  30  I.  =  0,23 
DO  20  N  = 1  , 200 

CALL  FADCK I  KAY ( (N-l )*NUM+l ) . NUM, IER) 
20  CONTINUE 


j  ' 

A 

l 

ft 

$ 


M 


TTT' 


STEr  THE  SCAN  I  VALVES  TO  THE  NEXT  POItT 

CALL  DOUTW ( I  DEV  1 CE , I STEP ) 

C  ALL  FDE1.Y  (  5 ) 

CALL  DOUTW ( 1  DEV  1 CE , 0 ) 

COMPUTE  THE  TOTAL  COUNTS  FOR  EACH  PORT  (OF  EACH  SCANIVALVE) 

AND  STORE  EACH  TOTAL  IN  THE  ARRAY,  Y. 

DO  30  M-NST ART, LIMIT 
Ml =M-NSTART+ 1 
DO  30  N =1.200 

Y(N1  ,L.)=Y(M1  , L)  +  IRAY(  (N-l  ) *NUM+  ( M-NSTART+ 1  ) ) 

30  CONTINUE 

RETURN  THE  SCAN  1 VALVES  TO  THE  HOME  POSITION 

CALL  DOUTW ( 1  DEV  ICE. 1 HOME) 

CALL  FDELY ( 5 ) 

CALL  DOUTW ( I  DEV  1 CE, 0  I 

DISPLAY  THE  AVERACE  SCANIVALVE  VALUES  (IN  COUNTS)  ON  THE  SCREEN 
NOTE:  THE  AVERACE  IS  COMPUTED  INSIDE  THE  WHITE  STATEMENT  ITSELF. 


DEC1N  TEMP.SV. 


TEMP . S V  >  >  >  BEC I N  CALPRES . SV . . 
CAl.PRES  .  SV  >>>  DEC  IN  NPARO.SV. 


END  NPARO.SV  >>>  CONTINUE  SCAN1.SV. 


; CONVERT  TO  PS I 


WRITE  FREE( 10) " 

CALL  FSWAP( “TEMP.SV “ ) 

WIU'IE  FREE(IO)"  END 

CALL  ESWAP (  "CAl.PRES .  SV "  ) 

WRITE  FREEdO)”  END 

CALL  ESWAP ( "NPARO.SV"  ) 

WRITE  FREER  10)" 

CALL  EOl'EN  (  3 ,  "TEMP"  ) 

READ  BINARY! 5)  TEMP 
CALL  ECLOSE ( 5 ) 

CALL  EOPEN  (  5,  "CAl.PRES  "  ) 

READ  Bl NARY ( 5 )  CAL 
CALL  ECLOSE ( 5) 

CALL  EOPEN ( 5 . “ATMPRES " ) 

READ  BINARY! 5)  ATM 
CALL  ECLOSE! 5) 

ATM=ATfR 3377/6U95 
WRITE  FREE! 10 ) "CAL" ,CAL 

TEHPR=TEHP+459.6  ; CONVERT  TO  RANK1NE 

TOTI’RES* (Y( 1 ,2)-Y( 1 .0) )*CAL^(Y( 1 , 1 )-Y< I .0) )+ATM 
STATPRES= (Y(4,2)-Y(4,0) ) .CAI./(Y(4. I ) -Y( 4 , 0 ) )+ATM 
rfL\Cll  =  SURT(ABS(3.*(  (  ( STATPRES/'TOTPRES ) ** ( -2 . /'7  .  )  )-l 
TRY I =5 . * ( ( ( STATPRES/TOTPRES ) ** ( -2 . /7 . ) )-l . ) 
TRY2=ABS(TRYI ) 

I F (TRY I  . EQ .TRY2 )  CO  TO  40 

WRITE  FREE! 10) "STATPRES  IS  CREATER  THAN  TOTPRES  ft" 
40  T1  =TEflPR/  (  1  +  (  FMACH**2  .  )  /5  .  ) 

A=SQRT! 1 . 4*1716. *T1 ) * . 3040 

FVEL=A*FMACII 

CALL  EOPEN ( 5 , "SCANI " ) 

WRITE  BINARY! 5)  TEMP.  ATM.  TOTPRES.  FMACH ,  FVEL 
CALL  FCLOSE! 5 ) 

C  REMOVE  THE  DEVICES 
C 

CALL  DIORM 
CALL  FADRM 

C 

STOP 

END 


)  )  ) 


-*»  m2* 


■y- ,  ■>  -/-V  -L-' v'vS-'v 

fr  - 1 1^ i.  lU-t. 


•  .V 
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*************************** 

*************************** 


xx  xx 
xx  “TEMP . FR“  xx 

XX  XX 
XX  MIKE  FOHEMAN  xx 
xx  xx 


*************************** 

*************************** 


"TEMP.FR"  COMPUTES  THE  CURRENT  ATMOSPHERIC 
TEMPERATURE  AT  THE  INLET  TO  THE  INDRAFT 
TUNNEL.  THIS  VALUE  IS  OUTPUT  TO  THE  DISPLAY 
SCREEN  AND  ALSO  WRITTEN  INTO  A  FILE  CALLED 
“TEMP" . 

DIMENSION  I  RAY ( 200 ) 


C 

0  THE  TEMPERATURE  PRODE  INTERFACES  WITH  CHANNEL  8 


NSTAKT  =  8 
LIMIT  =  8 

HUM  =  LI M1T-NSTART+ 1 


C 

C  SET  UP  THE  DEVICE 


0 

c 

0 


CALL  F A DST ( NSTART .LIMIT, 1ER) 
TAKE  200  SAMPLES 


DO  10  1=1. 200 

CALL  FADCT ( I RAY( I ) .HUM. IER) 
10  CONTINUE 


C 

•;  COMPUTE  THE  TOTAL  NUMDER  OF  COUNTS  FOR  THE  200  SAMPLES 
0 

TLllP-0  . 

DO  2<)  1  =  1, 200 
TEMP=TEMP+1RAY< I ) 

20  CONTINUE 
C 

C  CALCULATE  THE  AVERAGE  OF  THE  200  SAMPLES 
0 

TEHP= TEMP/200. 

C 

C  THE  TEMPERATURE  SIGNAL  ENTERS  THE  AMPLIFIER  IN  VOLTS. 

C  THE  AMPLIFIER  CIVES  THE  SIGNAL  A  OODB  GAIN  (FACTOR 

C  01  1000).  THEREFORE.  THE  SIGNAL  COES  TO  THE  A/D  CONVERTER  WITH  DIMENSIONS 

C  OF  MILLIVOLTS.  AFTER  CONVERSION.  THE  FACTOR  I63M.3  COUNTS/MI LLI VOLT 
C  IS  NECESSARY  TO  CONVERT  THE  COMPUTER  UNITS  (COUNTS)  BACK  TO  MILLIVOLTS: 


TEMP  =  TEMP/ 1 638 . 3  -.CONVERT  COUNTS  BACK  TO  MILLIVOLTS 
C 

C  THF.  "TEMP"  HAS  NOW  BEEN  CONVERTED  TO  MILLIVOLTS 
0  HOWEVER.  THE  TEMPEIIATUHE  EQUATION  REQUIRES  THAT 
C  THE  "TEMP"  BE  IN  MICROVOLTS: 


TEMP=TEMP* 1 000 . 


0  BECKMAN  COEFFICIENTS 
C 

CO  =  3  .  I  922787 F.+O  I 
Cl =3 . 5497020E-02 
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n  n  c>  n  nr>n  n  nn 


C2  =  -2 . 4700655E-07 


CALCULATE  TEMPERATURE 

TEMP=C0+C1*TEMP+C2*TEMP**2 
WRITE  TEMP  TO  SCREEN  AND  FILE  "TEMP" 
WRITE  FREE  (10) 

WRITE  FREE  (10)  "TEMPERATURE  (DECF)= 
CALL  FOPEN  (5."TEMP") 

WRITE  BINARY  (5)  TEMP 
CALL  FCLOSE  (5) 

DISCONNECT  THE  DEVICE 

CALL  FADRM 

STOP 

END 


* ,TEMP 
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*******************  ******** 
*************************** 


**  ** 

**  “CALPHES . FH”  ** 

**  ** 

**  MIKE  FOREMAN  ** 

**  ** 


*************************** 

*************************** 


“CALPRES . FR"  THIS  PROGRAM  COMPUTES  THE 
CALIBRATED  PRESSURE  TO  THE  SCAN I  VALVES. 

DIMENSION  I RAY ( 200 ) 

TIIE  CALIBRATION  PRESSURE  IS  CONNECTED  TO  CHANNEL  7. 

NSTAHT  =  7 
LIMIT  =  7 

NUM  -  L I  Ml T-NSTART+ 1 
SET  UP  THE  A/D  DEVICE: 

CALL  FADST(NSTART, LIMIT. IER) 

TAKE  200  SAMPLES  AT  CHANNEL  7 
DO  10  1=1, 200 

CALL  FADCT ( I  RAY ( I ) , NUM , I ER ) 

10  CONTINUE 

COMPUTE  THE  AVERAGE  OF  THE  200  SAMPLES 


Y  =  0. 

DO  20  1=1 .200 

Y  =  Y+FL0AT< IRAY(  1  ) ) 

20  CONTINUE 

Y=Y/20 O. 

CALIBRATION  PRESSURE  IN  INCHES  OF  HG  =  VOLTS* 10 
CAL=  (Y/1638.3)  M0. 

CALIBRATION  PRESSURE  IN  PS I  =  INCHES  OF  HC  *  (3377/6895) 


CAL=CAL# ( 3377 . /0895 . ) 

WRITE  FREE(IO)  CAL 

WRITE  THE  CALIBRATION  PRESSURE  TO  THE  SCREEN  AND  TO 
A  TEMPORARY  FILE  CALLED  “CALPRES”. 

WRITE  FREE(IO) 

WRITE  (10,25)  CAL 

25  FORMAT  ("CALIBRATION  PRESSURE:  *,F7.4,"  PSI " ) 

CALL  FOPEN  ( 5 . " CALPRES " ) 

WRITE  BINARY  (5)  CAL  ^ 

CALL  FCLOSE  (5) 

REMOVE  THE  DEVICE 


CALL  FADRM 


C  "PLOTP.FR"  IS  CALLED  BY  "SPOILER.FR"  TO  PLOT  SELECTED  VELOCITY 
C  PROFILES. 

C 

DIMENSION  X<25) .  Y(25) 

DIMENSION  1ST1 (3) , IST2I8) , IST3(8) , IST4IB) 

DIMENSION  I  STS ( 8 ) , 1ST6(9) , IST7I2I ) , IST8(3! ) 

C 

CALL  CURS  I Z( 4 ) 

DATA  I  ST 1/09, 47, 67/ 

DATA  I  ST 2/03 ,40,65.06,69,41.47,85/ 

DATA  1 ST3/06 , 40 , 65 , 86 ,69,41 .47,05/ 

DATA  I  S  I  4/05 ,40.02,77,03,41 ,47, 05/ 

DATA  1 ST5/G6 , 40 , 02 , 77 , 03 , 4 1 ,47,85/ 

DATA  1 ST6/B5 ,06.40,82.77,83,41,47, 05/ 

DATA  I  S  I  7/07 , 65 . 75 . 69 . 32 . 06 , 69 , 76 . 79 , 67 , 73 , 84 , 89 , 32 , 

*00 .02.79,70,73,76, 69/ 

DATA  I ST8/66 ,79.05,78.60.65,02,09.32,76,65,89,69,82,32, 

*06 ,69,76,79,67,73, 04 ,09 , 32 , 00 , 02 ,79,70,73,76, 69/ 

C 

CALL  FOPEN <  5 . “ PROF  I LE " ) 

READ  BINARY ( 5 )  NMAX ,  I  PROF ,  ISEL.  BICX 
READ  BINARY! 5)  (XU  )  ,  I  =  1  .  NMAX) 

READ  BINARY! 5)  ( Y (  I  )  ,  I  =  1  , NILAX ) 

CALL  FCLOSE ( 5 ) 

IFIBICX.LT . .5)  BIGX= . 5 
VR1TE  FREE!  1 0  ) Y(  I  »  , Y(NMAX)  ,BICX 
ACCEPT" ANY  NUI10LR  TO  CON  TINUE:  “  ,  IRAND 
C 

CALL  IN  ITT (960) 

CALL  B1N1TT 
CALL  CiiilS  1 Z  (  3  ) 

CALL  TERM (3, 1023) 

CALL  SLiriX!  100,900) 

CALL  SUHY(  150,650) 

CALL  XFRIK3) 

CALL  YFIUK5) 

CALL  XT I CS (10) 

CALL  YT 1 CS ( 10) 

CALL  DLiriX(O.O.BICX) 

CALL  DI. I I1Y  ( Y(  1  ) ,  Y ( NMAX) ) 

CALL  NPTS!  NILAX) 

CALL  CHECK (X.Y) 

CALL  DSPLAY(X.Y) 

CALL  CURS  I Z( 3) 

CALL  HOVABS! 10,525) 

CALL  VLABEL  (3, I ST1) 

CALL  MOVABS ( 300 , 725 ) 

IF! IPROF.EQ. 1 >  CALL  HLABEL( 2 1 , 1 ST7 ) 

I F ( IPR0F.EQ.2)  CALL  HLABEL(31 , IST8) 

CALL  MOVABS (420, 75) 

IFUSEL.EQ.I)  CALL  HLABEL(0,  IST2) 

I F ( ISEL.EO.2)  CALL  HLABEL( 8 , I ST3 ) 

IF(  ISEL.E0.3)  CALL  HLABEI.(8,  IST4) 

I  FI  ISEL. KQ. 4)  CALL  IILAR!.L(0,  IST5) 

IF!  ISEL.  ED.  5)  CALL  IILABEL ( 9 ,  I ST6  ) 

CALL  F I N ITT( © , 750 ) 

CALL  CIIRS I Z( 4  ) 

STOP 

END 
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c 


c 


c 


j******************************* 

******************************* 


** 

** 

** 

11 1 STOCRAM .  FR 

** 

** 

** 

** 

MIKE  FOREMAN 

** 

** 

** 

X*******X***XXX**XXX********XXX 

X*XX*X*X*XXXXX**XX**X*X***XX-kXX 

UISTOGHAM.FR  RECORDS  DATA  FROM  THE  LDV 
COUNTER  AND  PLOTS  A  HISTOCRAM  OF  THIS  DATA. 

REQUIRED  SUBROUTINES:  DI022  6FLI  BQ 

DIMENSION  1RAY< 1000) . IBIN( 16) ,FHEQ( 1000) 

DIMENSION  X<  101  )  .Y<  101  )  ,VF.l.(  1000) 

TYPE  "INPUT  SUUDEVICE  CODE  -  TWO  DIGITS  ONLY:" 

TYPE  "  CHANNEL  I  =  46" 

TYPE  "  CHANNEL  2  =  47" 

READ (11. 5)  I  SUB 
5  FORMAT (02) 

ACCEPT  "VELOCITY  CONSTANT?  “.CON 
ACCEPT  “MIXING  FREQUENCY?  " , FREQM 
10  ACCEPT  “NUMBER  OF  SAMPLES  (0  TO  STOP)  ?",NUM 
IF  (HUM. CT. 1000. OR.NUM.LT. 0)  GOTO  10 
IF  (Nim.EQ.O)  GO  TO  900 
CALL  DI01  I  (  I  SUB ,  I  RAY,  NUII,  I  ER) 

DO  100  J  =  1  , NUM 
DO  50  1=1,16 

I U!  N ( I )  =  I AND( IRAY( J) ,2**< 1-1  ) ) 

50  I F ( I B I W ( I > . N  E . 0 )  I B I N ( 1 )  =  1 
I  BN  =  0 

DO  60  1=0,9 
IBN=IBN+IBIN( 1+1 )*2**l 
60  CONTINUE 

I EXP  =  0 

DO  70  1=10,13 

I EXP= I EXP+ I B I N ( I + 1 )*2**( I -10) 

70  CONTINUE 

FREQ( J ) =32000 . /FLOAT < I BN*2** I  EXP ) 

VEL( J ) =CON * ( FREQ ( J ) -40 . +FREGH ) 

100  CONTINUE 

AVER=0 . 

A VE=0 . 

DO  I  1 0  1  =  1, NUM 
AVER  =  A VER+FREQ (  I  ) 

AVE= AVE+VEL( 1 ) 

110  CONTINUE 

A VER= AVER/HUM 
AVE= AVE/NUM 

WRITE  FHEE(IO)  "AVERAGE  DOPPLER  FREQUENCY=  ".AVER, “MHZ 
WRITE  FREK(IO)  "AVERAGE  VELOCITY*  ",AVE,“  M/S" 

ACCEPT  "DO  YOU  WISH  A  HISTOGRAM?  (YES=1,  N0=0)",AA 
IKAA.EQ.O)  CO  TO  999 
P1M=  I  FIX(AVE) 

J.I  =  fiM-5 

IF(.M.LT.O)  JJ  =  0 
KK  =  .I.I  + 1 0 
LL=  (.1.1*10)-! 


DO  200  1=1,101 
K= 1 +LL 

X! I ) = ( FLOAT ( K ) >/10. 

Y ( 1 ) =0. 

DO  150  J  = 1 , NUM 
1  I = I F1X( 10.*VEL( J ) ) 

I F < 1 1 .EO.K)  Y ( 1 )=Y( I )  +  l 
150  CONTINUE 

WRITE  FREE  (10)  X( I ) ,Y< I ) 

I W= I W+ I FIX( Y( I ) > 

I F  <  Y  < 1 ) . CT . Z )  Z  =  Y< I ) 

200  CONTINUE 

Z= FLOAT! ( I F IX! Z/ 1 00 . )*100)+I00> 

ik(z.ct. looo. )  z=iooo. 

AD  =  FL()AT(JJ> 

AC= FLOAT (KK) 

WRITE  FREE  (10)  AB.AC.Z 

CALL  FOPEN  (5.-HIST0" ) 

WRITE  BINARY  (5)  AD.AC.Z 
WRITE  BINARY  (5)  (X( I ) ,  1=1,101) 

WRITE  BINARY  (5)  (Y(I),  1=1,101) 

CALL  FCLOSE  (5) 

CALL  FSWAP  ( "HISTOCRAPH .SV“ ) 

WRITE  FREE(IO) 

WRITE  FREE( 10) 

WRITE  FREE! 10) 

WRITE  FREE! 10) 

WRITE  FREE  (10)  "  TOTAL  SAMPLES  PLOTTED 

CALL  DIORM 
999  STOP 
END 


<y 

fU 


“H1ST0CRAPH.FR"  USED  BY  H I STOCRAM . FR  TO  PLOT 
A  HISTOGRAM  OF  LASER  DATA. 

DIMENSION  X( 101 ) ,Y( 1 0 1 ) 

DIMENSION  I  ST 1(7),  IST2(21),  IST31I5) 

DATA  1  ST 1 /03 , 65 . 77 . 00 . 76 . 69 . 03/ 

DATA  I ST2/72 .73,03.04.79.71.02.65.77. 

*32.79,70,32,76. 60 , 06 , 32 . 60 , 65 , U4 ,65/ 

DATA  1 ST3/70 ,02,69,01,05,69,70,67.09.32,40,77,72,90,41/ 


CALL  FOPEN  (  5  ,  "11 1 STO"  ) 
READ  BINARY  <3)  AB.AC.Z 


READ  BINARY  (5)  (X( l ) . 
READ  BINARY  (5)  <Y< I) 
CALL  FCLGSE  (5) 


1=1,101) 

1=1,101) 


CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

SIOP 

EM) 


I N ITT( 960 ) 

BINITE 

TERM (3,1 023 ) 

V I1ARST  <  15.6,3) 
SI,  I  MX  (  100,900) 
SLII1Y<  150,650) 
Xl'RM  ( 2 ) 

YI'RiT  <  5 ) 

XT 1 CS ( 10) 

YIICSC 10) 

XfiTCS  (10) 

YMTCS (2 ) 
I)LIMX(AB.AC) 

Dl.  I  MY (0 . 0 ,  Z ) 
NI'TS(  101  ) 
CIILCK(X.Y) 
DSPLAY ( X, Y) 
M()VABS(  10,525) 
VLABEL(7 . IST1 ) 
MOVABS (350,725) 
IILAREL( 2 1  ,  I ST2 ) 
MOVABS (420. 75) 
1ILABF.L(  15,  I ST3  ) 
MOVABS  (  90') ,  130) 
DRV. AI1S  (900, 650) 
CIIRSIZ(4) 

FIN! TT(0,750) 


rn 
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